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ABSTRACT 

The spectral and temporal properties of the nonthermal emission of the nearby XRF 060218 in the 0.3-150 
keV band are studied. We show that both the spectral energy distribution and the light-curve properties suggest 
the same origin of the nonthermal emission detected by Swift BAT and XRT. This event has the longest pulse 
duration and spectral lag observed to date among the known GRBs. The pulse structure and its energy dependence 
are analogous to typical GRBs. By extrapolating the observed spectral lag to the Compton Gamma Ray Observatory 
0 CGRO ) B ATSE bands we find that the hypothesis that this event complies with the same luminosity-lag relation 
with bright GRBs cannot be ruled out at 2 a significance level. These intriguing facts, along with its compliance 
with the Amati relation, indicate that XRF 060218 shares the similar radiation physics as typical GRBs. 

Subject headings: gamma rays: bursts — methods: statistical 


1. INTRODUCTION 

X-ray flashes (XRFs), cosmic explosions with lower spectral 
peak energies (E p ) in vF v spectra than typical gamma-ray bursts 
(GRBs; Heise et al. 2001; Kippen et al. 2003), are thought to 
be the low-energy extension of typical GRBs (Lamb et al. 2005; 
Sakamoto et al. 2004, 2006; Cui et al. 2005). They may be 
GRB jets (uniform or structured) viewed at off-axis directions 
(e.g., Zhang et al. 2003, 2004; Yamazaki et al. 2004) or in- 
trinsically different events (Lamb et al. 2005; Soderberg et al. 

2005) . It is long speculated that long-duration GRBs originate 
from a relativistic jet emerging from a collapsing massive star 
progenitor (Woosley et al. 1993; Paczynski 1998; MacFadyen 
& Woosley 1999; Zhang et al. 2003), and associations of core- 
collapsing supemovae (SNe) with long GRB afterglows have 
been spectroscopically identified in a number of systems, in- 
cluding GRB 980425/SN 1998bw (Galama et al. 1998), GRB 
030329/SN2003 dhr (Stanek et al. 2003; Hjorth et al. 2003), 
and GRB 031203/SN 20031w (Malesani et al. 2004). The firm 
spectroscopic association of the nearby XRF 060218 (at z = 
0.0331; Mirabal et al. 2006) detected by Swift (Campana et al. 

2006) with the Type Ic SN 2006aj was established recently 
(Modjaz et al. 2006; Pian et al. 2006; Sollerman et al. 2006; 
Mirabal et al. 2006; Cobb et al. 2006; Soderberg et al. 2006). 
This suggests that the progenitors of both GRBs and XRFs are 
related to the death of massive stars, and XRFs are the low- 
energy extension of the more “standard” GRBs. 

Some empirical relations have been discovered from typical 
GRBs. It is interesting to verify whether XRF 060218 satisfies 
these relations as an ( approach to access its “standardness.” It 
has been reported that XRF 060218 complies with the isotropic 
energy versus the spectral peak energy (E iS0 ~E p ) relation derived 
from long GRBs (Amati et al. 2002) and other XRFs (Amati 
et al. 2006). The pulses in GRB light curves usually display 
a fast-rise-exponential-decay (FRED) shape, and it has been 
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found that the pulse width is related to the energy of the ob- 
servational band as oc E~ 0A (Fenimore et al. 1995; Norris et 
al. 2005). A correlation between the isotropic luminosity and 
the spectral lag (L^-t relation) of light curves was also dis- 
covered with typical GRBs as L iso oc r 118 (Norris et al. 2000). 
Both a)-E and L iso -r relations are related to the structure of light 
curves. Since most photons of an XRF are in the X-ray band, 
previous GRB missions did not observe the real light curves 
of XRFs. This makes it difficult to identify the temporal struc- 
ture and the spectral lag of XRFs. XRF 060218 is long and 
soft, and Swift Burst Alert Telescope (BAT) and X-Ray Tele- 
scope (XRT) simultaneously collected the data in the gamma- 
ray to X-ray bands. This makes it possible to measure its tem- 
poral structure and to examine whether it complies with the 
same o>~E and L iso -r relations derived from GRBs. In this Letter 
we focus on this issue. Throughout the Letter H 0 = 11 km 
s" 1 Mpc~\ Q m = 0.3, and 0 A = 0.7 are adopted. 

2. DATA 

XRF 060218 was detected with the Swift BAT on 2006 Feb- 
ruary 18.149 UT. It is a long burst, with a duration T 9Q ~ 
2000 s in the 15-150 keV band. Swift slewed autonomously to 
the burst, and the XRT and UV/Optical Telescope began col- 
lecting data 159 s after the burst trigger. 

XRF 060218 was an image trigger. The BAT event data 
lasted only until t ~ 300 s after the trigger. The survey data 
are used to derive the BAT light curve. XRF 060218 is very 
soft, with most of the emission in the BAT band being lower 
than 50 keV (Campana et al. 2006). In order to obtain a high 
level of signal-to-noise ratio of the BAT light curve, we use 
the light curve in the whole BAT band, i.e., 15-150 keV. The 
first orbit of the XRT data is fully in the Windowed Timing 
mode. We extract the light curves and spectrum of the XRT 
data with the Xselect package. The spectrum is grouped with 
the tool grppha, and the spectral fitting is carried out with 
the Xspec package. Following Campana et al. (2006) we fit 
the XRT spectrum in the first orbit with a model combining a 
blackbody component with temperature kT and a cutoff power 
law (F oc E~ r e~ El£c ) component. The absorption in both Milky 
Way 7 and the GRB host galaxy (A^ ost ) is incorporated. We 

7 The N h of our Galaxy (Dickey & Lockman 1990) for this burst is 
~1.1 x 10 21 cm" 2 . 
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Fig. 1. — (a) Unabsorbed light curves of the nonthermal gamma rays/X-rays in the energy bands of 15-150, 5-10, 2-5, and 0.3-2 keV, respectively. The fitting 
curves with eq. (1) are plotted. (b) Normalized light curves from the empirical model fitting. 


obtain /C 51 = 0.63 It o'o 3 x 10 22 cm" 2 , JcT = 0.122:°$* keV, 
T = 1.78^08, and E c — 2L8^ 5 7 8 keV (corresponding to 
E p ~ 5 keV), with reduced x 2 = 1.55 for 769 degrees of free- 
dom. We derive the unabsorbed light curves in the 0.3-2, 2- 
5, and 5-10 keV bands from the XRT data. As reported by 
Campana et al. (2006), a thermal component exists and is likely 
of different origin (e.g., the shock breakout emission associated 
with SN 2006aj; cf. Li 2006; Ghisellini et al. 2006b) from the 
nonthermal component. We therefore subtract the contribution 
of this component from the observed light curves. The AT and 
the radiation radii (R m ) evolve with time during the first orbit. 
We read the values of JcT and R BB from Campana et al.(2006) 
and calculate the light curves of this component in 0.3-2, 2- 
5, and 5-10 keV bands. Since the temperature of the blackbody 
component is below 0.2 keV, the light curves in the energy 
band higher than 2 keV is essentially not contaminated by the 
thermal component. In the 0.3-2 keV band, the derived light 
curve of the thermal component continuously increase with 
time, which could be well fitted by logT 03 _ 2 = (“10.69 ± 
0.09) + (0.66 ± 0.03) log t in cgs units. We thus subtract the 
contribution of this component from the 0.3-2 keV band light 
curve according to this fitting result. The derived nonthermal 
light curves in the three XRT bands as well as the one in the 


BAT band are shown in Figure 1 a. To clearly view the pulse 
width and spectral lag dependences with energy, we also plot 
in Figure lb the normalized light curves. The characteristics 
of the light curves along with the average photon energy (£) 
of each energy band are reported in Tables 1 and 2. 

3. RESULTS 

3.1. Spectral Energy Distribution 

One important question is whether the XRT nonthermal emis- 
sion and the BAT emission are of the same origin. In order to 
clarify this, we first study the spectral energy distribution (SED) 
of the burst using joint BAT-XRT data. Since this event is image 
trigger, the BAT event file contains only the data in the first 300 
s since the BAT trigger, XRT began collecting data 159 s after 
the B AT trigger. We thus only obtain simultaneous observations 
of the two instruments from 160 to 300 s with the event files. The 
joint-fit SED in this period is shown in Figure 2, which is well 
fitted by the BB+CPL model (Campana et al. 2006) with the 
following parameters: N^ ost = 0.65 :™ 3 x 10 22 cm"" 2 , JcT = 
0.13!™; keV, T = 1.56:™*, and£ c = 122^° keV (correspond- 
ing to E p ~ 54 keV), with reduced x 2 — 0.97 for 327 degrees of 
freedom. The E p strongly evolves with time, from 54 keV at the 


TABLE 1 


Temporal Structures of the Light Curves 


Band 

(keV) 

Peak 

(s) 

£*) 

(S) 

Rising Time 

(s) 

Decaying Time 
(s) 

<P 

E 

(keV) 

(!) 15-150 

.. 405(25) 

889(244) 

311(28) 

578(185 ) 

0.54(0.18) 

36.9 

(2) 5-10 

.. 735(9) 

1278 (45) 

475(12) 

803(35) 

0.59(0.03) 

6.9 

(3) 2-5 

.. 919(7) 

1707(40) 

624 (8) 

1084(34) 

0.58(0.02) 

3.1 

(4) 0.3-2 

.. 1082(13) 

2625(125) 

794 (14) 

1831 (112) 

0.43(0.03) 

0.7 


Note.— A ll the errors are derived by simulations and in the 1 a significance level. 
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TABLE 2 

Spectral Lags of the Light Curves 


Bands 

A E 
(keV) 

T peak 

(S) 

r CCF 

(5) 

Bands 

AE 

(keV) 

Tpeak 

(s) 

OM2) 

30 

330(26) 

249(37) 

UM3) 

33.8 

514(26) 

UM4) 

36.2 

677(28) 

518(70) 

(2H3) 

3.8 

184(11) 

(2)-(4) 

6.19 

347(16) 

173(25) 

(3)-(4) 

2.4 

163(15) 


Note. —A ll the errors are derived by simulations and in the 1 a significance 
level. 


beginning down to ~5 keV at later times. This is consistent with 
that reported by Campana et al. (2006) and Ghisellini et al. (2006a). 
The model fitting results are shown in Figure 2. One can observe 
that the BAT component is a good extrapolation of the XRT 
nonthermal component. This result implies that the nonthermal 
emissions detected by XRT and BAT are of the same origin. 



3.2. Pulse Width and Energy Dependence 

As shown in Figure 1, the light curves in different energy 
bands can all be modeled by a single FRED pulse. Kocevski 
et al. (2003) developed an empirical expression to .fit a FRED- 
like pulse, which reads 




(r+m-Rr+rfytr+l)] 


‘(I) 


where t m is the time of the maximum flux (F m )> t Q is the offset 
time, and r and d are the rising and decaying power-law indices, 
respectively. We fit the light curves with equation (1) and then 
measure the pulse width, rising and decaying times at the 
FWHM of the fitting light curves, and the rising-to-decaying 
time ratio (<p). The errors of these quantities are derived from 
simulations by assuming a normal distribution of the errors of 
the fitting parameters. The reported errors are at 1 a confidence 
level. The results are tabulated in Table 1. We show co as a 
function of E in Figure 3 (left). Apparently the two quantities 
are correlated. A best fit yields o) oc F“ a31±0 - 03 . The <p parameter 
ranges from 0.43 to 0.59. It is found that XRF 060218 roughly 
satisfies the same co-F relation (Fenimore et al. 1995; Norris 
et al. 2005), and its <p values are also well consistent with that 
observed in typical GRBs (e.g., Norris et al. 1996; Liang et 
al. 2002), although it has a much longer pulse width than other 
single-pulse GRBs. These results imply that XRF 060218 may 
be an extension of GRBs to the extremely long and soft regime. 

3.3. Spectral Lag and Energy Dependence 

The light curves shown in Figure la display a significant 
spectral lag (r), with soft photons lagging behind the hard 
photons, as usually seen in long GRBs (Norris et al. 2000; Yi 
et al. 2006). We illustrate this lag behavior with the intensity- 
normalized light curves in Figure lb. The light curves peak at 
405 ± 25, 735 ± 9, 919 ± 7, and 1082 ± 13 s, respectively, 
in a sequence of high-energy band to low-energy band as shown 
in Figure 3 (right). The best fit to the correlation between the 
peak time (r peak ) and the average photon energy yields 

log tpeak = 0.04 ± 0.04) - (0.25 ± 0.05) log F. (2) 

A simple estimate of the lags between any pairs of the four 
light curves obtains r peak = 163 ~ 677 s, being consistent with 
that shown in Gehrels et al. (2007). We note XRF 060218 


Fig. 2.— -BAT- XRT joint spectral energy distribution from 160 to 300 s since 
the BAT trigger. The BB+CutoffPL fitting model is also shown. 

becomes the new record holder of the long-lag, wide-pulse 
GRBs. The previous record holder was GRB 971208, with 
r ~ 58 s and co = 395 s (Norris et al. 2005). 

We also calculate the lags with the cross-correlation function 
(CCF) method. The errors of lags are evaluated by simulations. 
The results are also reported in Table 2. The lag derived by 
the CCF method (r CCF ) is strongly correlated with r peak but is 
systematically lower 8 than r peak (Fig. 4 [left]). A best fit gives 
t-ccf = (-100 ± 17) + (0.91 ± 0.08)r peak . 

The L iso -r relation was discovered with six bright BATSE 
GRBs (Norris et al. 2000), and the spectral lag was defined by 
the light curves in the 25-50 and 100-300 keV bands. We 
investigate whether the lag behavior of XRF 060218 is con- 
sistent with the L is0 -r relation. Since XRF 060218 is a soft XRF 
and the emission in the 100-300 keV band is too weak to 
derive a light curve, we assume that t peak of the light curve in 
the 100-300 keV band follows the r peak -F relation (eq. [2]) and 

8 We here derive r CCF from the peak of the CCF without considering the 
side lobe contribution of the CCF. A fit to the CCF with a cube or quartic 
function gives a larger lag by considering the side lobe contribution (Norris 
et al. 2000). However, this method strongly depends on the artificially selected 
range of CCF for the fitting. Since the light curves are a smooth pulse and 
their lags are significantly larger than the time bin, the peaks of CCFs are 
robust to estimate the lags. 



Fig. 3.-— Pulse duration (left) and the peak time (right) as a function of the 
average photon energy of the nonthermal emission. The solid lines in both 
panels are the best fits. 
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Fig. 4.— Left: Comparison of the spectral lags derived from the peak times 
and from the CCF method. The solid line is the best fit. Right: Isotropic gamma- 
ray luminosity as a function of spectral lag. The spectral lags of typical GRBs 
and GRB 980425 are calculated with the light curves in the 25-50 and 100- 
300 keV bands observed by CGRO BATSE. The lag of GRB 031203 is cal- 
culated with the light curves in the 20-50 and 100-200 keV bands. The gray 
band and the two dashed lines mark the best fits at the 1 and 2 a confidence 
level, respectively, and the solid line is the regression line for the six typical 
GRBs presented in Norris et al. (2000). 


perform the extrapolation. With the extrapolated f peak we then 
estimate r peak for the light curves in the 25-50 keV (average 
energy 30 keV) and 100-300 keV (average energy 200 keV) 
bands. We obtain r peak = 177 ± 16 s. Since r CCF is more re- 
liable, we use the r peak -r CCF relation (Fig. 4 [left]) to derive 
r CCF = 61 ± 26 s. This lag is used in the L-t relation analysis. 
Using the peak fluxes in the BAT and XRT band, we estimate 
L iso ~ 1.2 x 10 47 ergs s"\ 

Figure 4 (right) shows the L iso -r relation derived by Norris et 
al. (2000) compared against XRF 060218 as well as two other 
nearby GRBs, 980425 and 031203. The data of the previous 
GRBs are taken from Norris et al. (2000) and Sazonov et al. 
(2004). 9 The gray band and the two dashed lines mark the best 
fits at the 1 and 2 a confidence level, respectively, and the solid 
line is the regression line for the six GRBs that were used to 
draw the L-t correlation; i.e., log L iso = (50.22 ± 0.32) - 
(1.21 ± 0.21) log r (errors arc at the 1 a level). We can see that 
XRF 060218 is definitely inside the 2 a region and is marginally 
at the 1 a region boundary. Therefore, the hypothesis that XRF 
060218 follows the L-t relation cannot be ruled out at the 2 cr 
significance level. We caution that the r is inferred from the 
extrapolation of the f peak -£ relation. This introduces uncertainties 
in deriving the lag. The other two nearby GRBs, 980425 and 
031203, are out of the 2 a region, which are identified as sig- 
nificant outliers of this relation (e.g., Sazonov et al. 2004). 

9 The peak photon fluxed in the 50-300 keV bands of these GRBs are 
converted to energy fluxes with the spectral index presented in Friedman & 
Bloom (2005), and their L^-values are recalculated with the cosmological 
parameters used in this Letter. 
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4. DISCUSSION AND CONCLUSIONS 

We have investigated the nonthermal emission of XRF 
060218. The early SED of this event from 0.3-150 keV ob- 
served by BAT and XRT suggests that the nonthermal emission 
detected by the two instruments are the same component. By 
subtracting the contribution of the thermal emission we derive 
the light curves of the nonthermal emission. They arc composed 
of a broad single pulse, and the energy dependences of the 
widths and the rising-to-deeaying-time ratio of the pulses are 
roughly consistent with those derived in typical GRBs. The 
light curves show significant spectral lags, with a well-defined 
peak time sequence from high-energy band to low-energy 
bands; i.e., oc L” 0-25 ” 0 - 05 . We infer the spectral lag in the 
BATSE bands and find that the hypothesis that this event com- 
plies with the L iso -r relation with typical GRBs cannot be ruled 
out at the 2 a significance level. 

These intriguing facts, along with its compliance with the 
Amati relation, strongly suggest that GRB 060218 is a standard 
burst at the very faint, long, and soft end of the GRB distri- 
bution. Since all these relations concern the temporal and spec- 
tral properties of emission, they are likely related to the radi- 
ation mechanisms. The results therefore imply that XRF 
060218 and other XRFs may share the similar radiation physics 
(e.g., synchrotron or inverse Compton scattering in internal 
shocks; Meszaros 2002, 2006; Zhang & Meszaros 2004; Piran 
2005) with harder GRBs. 

As discovered by Norris (2002), the proportion of long-lag 
bursts within long-duration bursts increases from negligible 
among bright BATSE bursts to -50% at the trigger threshold, 
and their peak fluxes are -2 orders of magnitude lower than 
those of the brightest bursts. This argues that they are intrinsically 
underluminous. Taken together with the fact that three nearby 
GRBs, 980425, 031203, and 060218, arc long-lagged and un- 
derluminous, an intuitive speculation is that long-lag bursts are 
probably relatively nearby (e.g., Norris et al. 2005). The local 
GRB rate of these GRBs should thus be much higher than that 
expected from the high-luminosity GRBs (Liang et al. 2006; see 
also Cobb et al. 2006; Pian et al. 2006; Soderberg et al. 2006). 
A possible scenario to explain their wide-pulse, long-lag, and 
underluminous features is die off-axis viewing angle effect (e.g., 
Nakamura 1999; Salmonson 2000; Ioka & Nakamura 2001). 
Another scenario is that these features are intrinsic, being due 
to their lower Lorentz factors (Kulkami et al. 1998; Woosley & 
MacFadyen 1999; Salmonson 2000; Dai et al. 2006; Wang et 
al. 2006). They might be from a unique GRB population (Liang 
et al. 2006), having a different type of central engine (e.g., neu- 
tron stars rather than black holes) from bright GRBs (e.g., Maz- 
zali et al. 2006; Soderberg et al. 2006). 
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